Dynamical quantum phase transitions 
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A sweep through a quantum phase transition by means of a time-dependent external parameter (e.g., pressure) 
entails non-equilibrium phenomena associated with a break-down of adiabaticity: At the critical point, the 
energy gap vanishes and the response time diverges (in the thermodynamic limit). Consequently, the external 
time-dependence inevitably drives the system out of equilibrium, i.e., away from the ground state, if we assume 
zero temperature initially. In this way, the initial quantum fluctuations can be drastically amplified and may 
become observable - especially for symmetry-breaking (restoring) transitions. By means of several examples, 
possible effects of these amplified quantum fluctuations are studied and universal features (such as freezing) are 
discussed. 



I. MOTIVATION 

In many cases, the properties of the ground state (or, at fi- 
nite temperatures, the thermal equilibrium state) are sufficient 
for understanding the behavior of a given condensed-matter 
system. In many other scenarios, however, the actual quantum 
state shows significant deviations from the ground state and 
non-equilibrium phenomena have to be taken into account. 
Even if the system starts near the (initial) ground state, it may 
depart during the evolution - for example, due to an external 
time-dependence. 

For example, let us consider an explicitly time-dependent 
Hamiltonian with a discrete (non-degenerate) spectrum 
H(t)\\j/ n (t)) = E„(t) | Yn(f)). After inserting this instanta- 
neous eigenvector expansion into the Schrodinger equation 
and some approximations, one obtains the adiabatic expan- 
sion 
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I.e., if the system started in its ground state \tyo{to)}, its ac- 
tual quantum state |y(/)) will remain near the instanta- 
neous ground state | y/b(0) if the non-adiabatic corrections are 
small, i.e., if (\ff„\H\\jfQ) <C (E„ — Eq) 2 . The latter condition 
compares the rapidity of the external time dependence H with 
the internal response times given by the energy level spacings 
AE„ = E n — Eq. If the response time of the system is short 
enough such that it can adapt to the external variation, it will 
stay near the ground state - if the externally imposed change 
is too fast, however, the system cannot respond and will de- 
part from the instantaneous ground state \\j/o(t)), leading to 
the creation of excitations \ \j/„>o), e.g., quasi-particles. 

This observation directly indicates the close ties between 
(quantum) phase transitions jj] and non-equilibrium phenom- 
ena: At the critical point, at least some of the energy levels 
converge AE„ = (E n — Eq) l and hence the associated re- 
sponse times diverge AEAt ~ h, which makes it very easy to 
depart from equilibrium, even for very slow external time de- 
pendences H, see, e.g., J2H5I Of course, the resulting non- 



equilibrium phenomena do not just depend on the energy lev- 
els AE„ = (E n — Eq), the matrix elements (y n \H |Vb) are a ls° 
very important. 



II. FIRST EXAMPLE 

Let us start by considering a very simple example: a homo- 
geneous two-component Bose-Einstein condensate |6]. In the 
dilute-gas limit, the Hamiltonian density reads 
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where = (¥1,^2) denote the field operators of the two 
components, m the mass of the condensed particles, and g a t, 
the coupling matrix in the ^-scattering approximation. Using 
the mean-field approximation and Madelung spilt *F a sif fl = 
v /p^exp(/4>„) with the densities p„ > and the phases 4>„, 
we obtain two decoupled Goldstone modes due to the bro- 
ken [7(1)017(1) symmetry 4>± = (<54>i ± <54> 2 )/\/2. The 
hard density mode <I> + corresponds to variations of the to- 
tal density p = Pi +P2 while the soft spin mode <t> is as- 
sociated to the density difference p\ — p 2 - The propaga- 
tion velocities of the two modes are given by the eigenval- 
ues g± = {g u +g 22 ) 2 /4±^(gu -g22) 2 / 4 + gi2 ofthe cou- 
pling matrix g a j,. Depending on the relation between the inter- 
component repulsion g\2 = g2i > and the intra-component 
couplings gn,g22 > 0, there is a phase transition when g_ 
changes its sign. If the repulsion within the same component 
gi i)g22 dominates, the condensate will be in the homogeneous 
mixed state p\ = In case the particles of different com- 
ponents repel each other more strongly, the inhomogeneous 
phase separated state (either p\ = or P2 = 0) will be ener- 
getically favorable, see Fig.[T] 

Now, let us study the sweep through this phase transition by 
means of a time-dependent g 12(f) starting in the mixed state. 
The hard density modes <t> + will basically remain unaffected 
by the transition, but the spin modes <t>_ become arbitrarily 
soft when approaching the critical point g_ = and finally 
become unstable on the other side. Hence these modes are 
the most interesting ones. In such a sweep, there will be four 
major stages: 



2 




FIG. 1: Sketch of the level structure from Eq. (O near the critical 
point g_ = separating the mixed condensate from the phase sepa- 
rated state. The two solid lines denote these two competing ground 
states which cross at g_ = and the quasi-particle excitations are 
depicted by dashed lines. 

• cooling Since approaching the critical point implies a 
reduction of the quasi-particle energies AE„ = E„ — Eq 
of the spin modes, their temperature (if non-zero ini- 
tially) will drop as long as we are still approximately in 
equilibrium (i.e., adiabatic). 

• freezing Close enough to the critical point, however, 
the spin modes become too soft, i.e., slow, and thus can- 
not adopt to the externally imposed variation g 12(f) any- 
more, i.e., they freeze. 

• squeezing Even if we started in the initial ground 
state, this freezing of the modes implies a squeezing 
(and thus amplification) of the quantum fluctuations due 
to the breakdown of adiabaticity, cf. Eq. (Q} and Fig. [2] 

• re-heating Finally, far enough on the other side of the 
phase transition, the system starts to "realize" that it 
is not in its true ground state and the spin modes star 
to grow exponentially - with the frozen and amplified 
quantum fluctuations as initial seeds. 

Interestingly, this sequence: cooling — > freezing — > squeezing 
— > re-heating exhibits striking similarities to cosmic inflation. 
Indeed, one might speculate whether cosmic inflation was not 
a real expansion of the universe over many orders of mag- 
nitude, but instead represent our distorted view of an early 
cosmic phase transition [0. 



III. HORIZON ANALOGUE 

Let us study the general behavior mentioned in the previ- 
ous Section in some more detail. To this end, we have to con- 
sider the equation of motion for the spin modes. Of course, 
we could derive it directly from the original Hamiltonian (f2]). 
However, for the sake of universality, let us base our discus- 
sion on general arguments instead. Assuming that the spin 
modes can be described by a linearized low-energy effective 




FIG. 2: Schematic of the squeezing (amplification) of the quantum 
fluctuations. Each (linearized) quasi-particle mode corresponds to a 
harmonic oscillator, whose potential is depicted as a dashed curve. 
Initially (left picture), the mode is in its ground state (solid curve). 
Approaching the critical point, the modes become softer and hence 
the potential gets more shallow (right picture). As long as this hap- 
pens slow enough [cf. Eq. (Q}], the quantum state evolves adiabat- 
ically and stays near the instantaneous ground state (dotted curve). 
After freezing, however, the modes cannot adapt to this change any- 
more and their quantum state (solid curve in right picture) deviates 
from the ground state (dotted curve). In view of the linearity of the 
equations of motion, the quantum state evolves into a squeezed state 
(e.g., the wave function is still Gaussian). This squeezing effect be- 
comes even more pronounced if the quadratic potential turns over, as 
in Fig.|9] The direction in which the quantum fluctuations will be am- 
plified depends on the system under consideration. For a symmetry- 
breaking transition (where the potential turns over, cf. Fig. [3), for 
example, the quantum fluctuations typically grow in the direction of 
symmetry-breaking. 

action, we write down the most general form for scalar Gold- 
stone modes in a homogeneous and isotropic medium |0] 

where the two remaining factors a and j3 are functions of the 
external parameter g_ (t) and have to be determined from the 
underlying system. Accordingly, the effective energy (Hamil- 
tonian) density of these Goldstone modes reads 

^eff = ^(afe(0]n 2 + )3fe(r)](V4>) 2 ) . (4) 

Before the transition, the spin modes are stable and hence a 
and j3 must be positive. After the transition, however, the spin 
modes become unstable and thus at least one of the two pa- 
rameters has to change its sign, i.e., it must vanish at the criti- 
cal point. As a result, the propagation velocity c 2 = aj3 must 
also go to zero when approaching the transition, which pre- 
cisely reflects the fact that the spin modes become infinitely 
soft. 

The consequence of such an evanescent propagation ve- 
locity c(f) can nicely be explained by the analogy to grav- 
ity/cosmology: For a constant speed dc/dt = 0, quasi-particle 
excitations may propagate arbitrarily far through the sample 
(in the absence of damping etc.) if we wait long enough. If 
c(t) decreases sufficiently fast, however, the excitations may 
not travel infinitely far - even if one waits long enough and 
takes an infinite time to reach the critical point - but merely a 
finite distance, which is an exact analogue of a cosmic horizon 

00 

Ar(f) = J dt'c{t'). (5) 
t 
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The horizon size Ar(f ) corresponds to the maximum distance 
covered by quasi-particles starting at time t and hence mea- 
sures the range of causal connection [IS]. Of course, in a real 
experiment, the critical point is crossed after a finite time, 
which should then be used as the upper limit (with the physics 
conclusions remaining unchanged). 

Two points at distances larger than the horizon Ar(t) can 
no longer exchange information or energy, i.e., the horizon 
describes a dynamical loss of causal connection. Since main- 
taining equilibrium requires causal connections between all 
points of the sample (e.g., in order to equilibrate local density 
variations), this loss of causal connection indicates the depar- 
ture from equilibrium and the breakdown of adiabaticity. As 
mentioned before, this typically entails an amplification of the 
quantum fluctuations, which can be understood in the follow- 
ing way: Since the analogue horizon shrinks steadily 

|mo=|/^v(o=-c ( o ! (6) 

t 

it will engulf every wavelength X after a given time. There- 
fore, all modes X will pass through the three main stages of 
their evolution: 

• oscillation Initially X <C Ar{t), the modes do not no- 
tice the horizon and oscillate almost freely (though with 
a decreasing frequency, cf. the cooling in the previous 
Section). 

• horizon crossing At some point of time, the horizon 
closes in X ~ Ar(t ) and inhibits further oscillations be- 
cause crest and trough of a wave are separated by the 
horizon and hence cannot exchange energy anymore. 

• freezing & squeezing After that, the modes are 
frozen and cannot respond to the external change any- 
more, i.e., the initial Gaussian ground state will turn 
into a squeezed state representing the amplification of 
the quantum fluctuations, cf. Fig. [2] 

According to our standard model of cosmology, precisely the 
same mechanism occurred in the early universe (during in- 
flation) and is responsible for the generation of the seeds for 
structure formation out of the initial quantum vacuum fluctua- 
tions. Traces of these amplified quantum vacuum fluctuations 
can still be observed today in the anisotropies of the cosmic 
microwave background radiation. 

IV. UNIVERSALITY? 

Now, after having studied this specific example, one is 
lead to the question of how universal these phenomena re- 
ally are. In order to study this point, let us have a look at 
the effective energy landscape in terms of the order parame- 
ter x = Pi/(pi + P2). It is restricted according to < x < 1 
with the boundary values representing the phase separated re- 
gion while the mixed state is given by x = 1/2. Inserting this 
ansatz into the Hamiltonian (O, the relevant part of the energy 




FIG. 3: Sketch of the energy landscape associated to Eq. The 
filled black dot denotes the initial ground state (mixed condensate) 
and the hollow dots correspond to the final ground states (phase sep- 
aration). In view of the jump in the ground state structure, the phase 
transition is of first order - though it is a somewhat untypical exam- 
ple, since the energy landscape becomes flat (in the relevant direc- 
tion) at the critical point (middle panel). 




FIG. 4: Sketch of the energy landscape associated to a typical first- 
order phase transition. Again the filled/hollow black dot denotes the 
initial/final ground state. In this case, the vicinity of the initial ground 
state remains (locally) stable at the critical point (second panel from 
the left) and the linearized quasi-particle modes do not become ar- 
bitrarily soft at this stage. Typically, this happens much later, at the 
point of instability of the initial ground state (right panel). 



functional reads g_ (x — 1/2) 2 , which is plotted in Fig. [3] Even 
though the phase transition is of first order (on the mean-field 
level), it is a somewhat untypical example, because the energy 
landscape becomes flat at the critical point. Therefore, adding 
a small non-linearity such as (x — 1 /2) 4 would turn it into a 
transition of second order. 

In contrast, the energy landscape of a typical first-order 
transition is sketched in Fig. [4] In this case, the vicinity of the 
initial ground state remains stable even after the transition and 
there is no local linearized mode connecting the two compet- 
ing vacua before one reaches the point of instability, cf. Fig. [4] 
right panel. Therefore, the previous arguments (softening of 
modes etc.) do not apply to the critical point, but they could 
be applied to the point of instability (which occurs later). 

The aforementioned instability of the transition in Fig. [3] 
towards the transformation into a second-order transition al- 
ready suggests that the phenomena discussed in the previous 
Sections could also be relevant for transitions of higher order. 
Fig. [5] depicts the energy landscape of a typical symmetry- 
breaking second-order phase transition. As suggested by this 
picture, there will be some modes which become arbitrarily 
soft at the critical point and unstable afterwards, i.e., during 
the sweep through such a transition, one generally reproduces 
the sequence: cooling — s- freezing — > squeezing — > re-heating. 
Transversing the critical point in the opposite direction (i.e., 
a symmetry-restoring transition), there are still modes which 
become arbitrarily soft - but they do not necessarily become 
unstable. Consequently, the shorter sequence: cooling — ^ 
freezing — > squeezing does also apply here, but the re-heating 
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FIG. 5: Sketch of the energy landscape associated to a typical 
symmetry-breaking second-order transition. At the critical point 
(middle panel), the local curvature in the vicinity of the initial 
ground state (filled dot) changes its sign and hence the associated 
quasi-particle modes become arbitrarily soft at the transition and 
turn unstable afterwards. The time-reversed sweep corresponds to a 
symmetry-restoring transition. In this case, the ground state (hollow 
dot) does not become unstable, but still the associated quasi-particle 
modes become arbitrarily soft at the transition. 




FIG. 6: Typical evolution of the dispersion relation where the insta- 
bility occurs at a finite wavenumber k CIIt . Before the transition (solid 
curve), the frequencies are non-negative. At the critical point (dashed 
curve), the "roton" dip touches the axis and afterwards, the modes in 
the vicinity of fc cr ; t start to grow (dotted curve). 
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FIG. 7: Typical evolution of the dispersion relation where the insta- 
bility occurs at k = since the quasi-particle gap - i.e., their effective 
mass m 2 in Eq. ((7} - changes sign. 
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FIG. 8: Typical evolution of the dispersion relation where the in- 
stability occurs at k = since the quasi-particle stiffness c 2 of the 
Goldstone modes changes sign, cf. (3]l. 



part requires further considerations. 

Of course, having established that some modes become ar- 
bitrarily soft when approaching the critical point does not 
mean that they behave like the spin modes in Eq. (0). For 
the example in Sec. HH the typical length scale of the am- 
plified quantum fluctuations was set by the horizon size and 
thus the sweep rate dc/dt. This is not always the case. In 
order to explore the different possibilities, let us study quasi- 
particle modes which become arbitrarily soft when approach- 
ing the critical point (and unstable afterwards). Their rele- 
vant features should be encoded in the quasi-particle disper- 
sion relation co(k). Assuming invariance under time-reversal 
and reflection as well as (spatial) rotations lErjll . it simplifies to 
(0 2 (k 2 ). The critical point then occurs when (0 2 (k 2 ) starts to 
dive below the vertical axis. This could happen at non-zero k, 
see Fig. [6] or at zero k. In the latter case, we may employ a 
Taylor (i.e., low-energy) expansion similar to Eq. <£3j 



eo 2 Ck 2 )=m 2 +c 2 k 2 + 0{k 4 ). 



(7) 



Thus there are two options for an instability: either m 2 
changes its sign, see Fig. [71 or c 2 changes its sign while m = 0, 
see Fig. [8] Note that the case of c 2 changing its sign in the 
presence of a positive m 2 would bring us back to the case in 
Fig. [6] 

An example for the scenario in Fig.|6]is the extended Bose- 
Hubbard model (cf. Sec. IVIt where the dispersion relation 
may possess a "roton" dip which dives below the vertical axis 
at the superfluid-supersolid phase transition. In this case, the 



typical length scale of the amplified quantum fluctuations is 
naturally set by the characteristic wavenumber k cr i t at the "ro- 
ton" dip, which determines the period of the supersolid order 
and is independent of the sweep rate etc. However, spatial 
modulations around this wavenumber A: cr i t (which could in- 
duce defects in the supersolid crystal, for example) may de- 
pend on the sweep rate. 

The instability at small k sketched in Fig. |7j occurs in the 
example discussed in Sec. IVIII In contrast to Fig. [6] the 
characteristic length scale cannot be read off the critical dis- 
persion relation (dashed curve) in general. Instead, it is set 
by the sweep rate and/or the final state 1271, d epending on 
the system parameters. For the experiment 111311 discussed in 
Sec. IVIII the final state determines the characteristic length 
scale. In other regions of parameter space (e.g., with a much 
slower quench), however, the sweep rate may yield the dom- 
inant contribution. In this situation, a quantum version of 
the Kibble-Zurek ifTii [l5ll scaling arguments j 1(1 may be ap- 
plied: The time-dependent correlation length |(f) <*; \/m(t) 
is cut off at the time ? when adiabaticity breaks down, i.e., 
? ~ Response (?) x 1/A£(f). This saturated correlation length 
B, = B, (?) then determines the typical distance of the topologi- 
cal defects. 

Finally, the case of Goldstone modes (which become unsta- 
ble upon crossing the transition) is sketched in Fig. [8] For ex- 
ample, the system discussed in Sec.[n]falls into this category. 
Note that Fig. [3] describes the behavior of that system under 
large but homogeneous (i.e., k = 0) displacements, whereas 
the dispersion relation in Fig. [8] corresponds to small (lin- 
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earized) inhomogeneous displacements with various k. In this 
case, the concept of a horizon analogue sketched in Sec. [Til] 
can be applied. In the vicinity of the critical point (dashed 
line in Fig. [8]), the time-dependent horizon size in Eq. (0 de- 
termines the spectrum j7] and the typical length scales of the 
amplified quantum fluctuations, but after the transition (dot- 
ted line in Fig. [B), the final state typically induces an addi- 
tional length scale corresponding to the fastest growth of the 
modes (negative minima of dispersion curve). Note that the 
Kibble-Zurek scaling arguments sketched above do not apply 
in this situation: Modes with different k cross the horizon in 
Eq. Q and thus freeze (i.e., become non-adiabatic) at differ- 
ent times - hence, there is no unique time f and consequently 
no saturated correlation length t, = ^ (t ). 

V. BOSE-HUBBARD MODEL 

Now, after having discussed some general aspects, let us 
turn to a few explicit examples. A prototypical example |01 for 
a second-order phase transition is the Bose-Hubbard model 

H = J(t)Y,M a paia p + y £(4) 2 a 2 ■ (8) 
a/3 a 

As usual, da, dp are the (bosonic) creation/annihilation oper- 
ators for the lattice sites a and j3, respectively. Hence, the 
first term describes hopping between the lattice sites a and 
j3, where the lattice structure (e.g., cubic) is encoded in the 
matrix M a p and the tunneling rate 7(f) is chosen to be time- 
dependent. The second term describes the fact that two bosons 
repel each other if they are on the same lattice site, where 
U denotes the associated energy penalty. For integer filling 
n = (d a d a ) = («a) € N, the system undergoes a second-order 
quantum phase transition l28ll at a critical value of the hopping 
rate 7 cr ; t = &(IJ jn). For J ^> U/n, the first term dominates 
and the bosons can move freely on the lattice, which corre- 
sponds to the superfluid phase 

/>E//n^|¥) sf oc (£4) l°>> (9) 

where N = £ a n is the total number of bosons on the lattice. 
In the opposite limit 7 < U/n, the second term dominates and 
causes the particles to lock in the lattice sites, i.e., they are not 
free to move across the lattice anymore. Hence this phase is 
called the Mott insulator state 

J « U/n ^ |¥> Mott « 0(4)" |0) . (10) 

a 

In the Mott insulator phase, the ground state is separated by a 
finite gap from the other excited states. The superfluid state, 
on the other hand, is associated to a certain phase and thus 
breaks the 17(1) symmetry d a —> e'^d a of the Hamiltonian 
(0. As a result, this state supports gap-less Goldstone modes, 
which are the lattice phonons. 

Now, let us consider the following scenario: starting deep 
in the superfluid phase J ^> U /«, we decrease the hopping 
rate J(t) and thereby sweep through the superfluid-Mott phase 



transition far into the Mott regime / <[/ /«. If we did this in- 
finitely slowly, we would end up in the Mott state dTOb . A 
very rapid sweep, on the other hand, would leave the system 
no time to respond and thus it would stay in the initial coher- 
ent state (|9]l. For a finite sweep velocity dJ/dt, one would 
expect to end up somewhere in between the two states. A 
good marker for distinguishing the two regimes is the (final) 
number variance A 2 (n a ) = (h a ) — (««) 2 - The coherent state 
generates Poissonian number statistics and these fluctuations 
are large A 2 (n a ) = n in the superfluid state J^$> U jn, whereas 
they vanish deep in the Mott state ( fTob . 

Because the Hamiltonian © cannot be diagonalized ana- 
lytically, some approximations are necessary for deriving ex- 
plicit results J^,[l3]. Here, we assume a large filling n> 1, 
which facilitates an expansion into (inverse) powers of n via 
n a =n + 8h a + €?(n°) where the fluctuations 8n a = ff{^Jn) 
can be treated analytically. These fluctuations correspond 
to the Goldstone modes mentioned earlier and their propa- 
gation speed (in the continuum limit) can be estimated via 
c 2 = £ 2 JUn, where I is the lattice spacing. Thus, if J de- 
creases fast enough, we get a horizon analogue as discussed 
in Sec. mil 

Let us study a few examples for the temporal behavior of 
the sweep: For an exponentially decaying hopping rate /(/) = 
/oexp{ — yt}, a horizon emerges for all values of y and hence 
the Goldstone modes freeze, cf. Sec. HH] The frozen number 
fluctuations are given by 

1 — p~ 2liv 

A 2 (n a ) = {hi) (n a ) 2 = (8n 2 a ) = n ^ , (1 1) 

where the adiabaticity parameter V = Un/y = n/yis given by 
the ratio of the internal time scale (chemical potential fi=Un) 
over the external sweep rate y and thus measures the rapidity 
of the sweep. In agreement with the previous arguments, we 
recover the limiting cases A 2 (n a ) f n for V 4 and A 2 {n a ) \. 
for V t °°- If the hopping rate merely exhibits a slow poly- 
nomial decay J(t) oc t~ x with x < 2, no horizon emerges and 
the number fluctuations do not freeze at a finite time but oscil- 
late forever. The boundary case J(t) °< t~ 2 is quite interesting: 
Since it marks the limit for horizon formation, the modes show 
a behavior in between freezing at a finite value and eternal os- 
cillation, i.e., they decay down to zero quite rapidly. This situ- 
ation is analogous to a damped harmonic oscillator, where the 
critical damping [corresponding to 7(f) °= f ~ 2 ] marks the bor- 
der between the over-damped (horizon formation and freez- 
ing) and the under-damped case (eternal oscillation). Hence 
this sweep dynamics J(t) ^ t~ 2 would be the optimal choice 
for approaching the Mott state with A 2 (n a ) \. most quickly 



VI. EXTENDED BOSE-HUBBARD MODEL 

If we extend the Hamiltonian (O and include inter-site in- 
teractions given by V a p 

H = Y J { T a^ a dp+V a ph a np) , (12) 
ap 
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the phase diagram becomes much richer and now includes su- 
persolid states etc. The superfluid-supersolid phase transition 
is a nice example for the scenario in Fig. [6] where the disper- 
sion relation develops a "roton" dip which touches the ver- 
tical axis at the critical point (at zero temperature, cf. ITlll ). 
This implies an instability towards density modulations at a 
critical wavenumber k CI it, which corresponds to the period of 
the supersolid order parameter. As before, the modes be- 
come arbitrarily soft when approaching the transition - which 
entails interesting non-equilibrium effects. For example, in 
the presence of an external phase gradient (i.e., condensate 
flow), the branch of the "roton" dip corresponding to quasi- 
particles propagating against the condensate flow is closer to 
the axis, i.e., softer, than the other branch. Hence, approach- 
ing the transition too fast will create non-adiabatic excitations, 
which are stronger for these (softer) modes, because their gap 
is smaller, cf. Eq. (fl}. As a result, non-equilibrium phenom- 
ena may induce a current (stemming from the amplified quan- 
tum fluctuations) in the opposite direction - which might even 
compensate the condensate flow fllll . 



VII. SPINOR (5=1) CONDENSATE 

As an example for the behavior sketched in Figs. [5] and [7] 
let us consider a three-component Bose-Einstein condensate, 
which can be represented by an effective spinor field operator 
*£ = (^ Xy W y ,4' z ) and the Hamiltonian density 13 



— \ 2m 2~ — ~ 



4in>?c 



<7out <9 a 



(13) 



The first term is the same as in the case of one or two com- 
ponents (O and contains the pseudo-spin independent kinetic 
term V / (2m) and coupling cq. The second term C2 corre- 
sponds to spin-1 scattering channel and favors the ferromag- 
netic phase with maximum magnetization F_ = x ^ 
due to C2 < 0. Finally, the third term q > denotes the 
quadratic Zeeman shift and would be minimized by placing 
all particles in the z-component, i.e., it favors the paramag- 
netic state ()£ X = 0. The (second-order) phase transi- 
tion separating the two regimes occurs at the critical point 
<7crit = 2\c 2 \p. 

Motivated by a recent experiment [13], we study the sweep 
from the paramagnetic to the ferromagnetic phase, cf. Fig. [9] 
Initially, the system stays in the potential minimum at zero 
magnetization, but after the quench, the initial minimum at 
F_ = turns to a local maximum and the system starts to roll 
down the potential hill, cf. Fig. [9] Since this is a symmetry- 
breaking phase transition, the direction of descent can be ar- 
bitrary and is set by the initial quantum (or thermal) fluctua- 
tions. Thus the spontaneously formed (measurable) directions 
of the magnetization F_ provide a direct indicator for the ini- 
tial quantum fluctuations, which were strongly amplified in 
that process. One result of these amplified quantum fluctua- 
tions could be the creation of topological defects (which cor- 
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FIG. 9: Sketch of the energy landscape associated to Eq. d 1 3 b - Ini- 
tially, the system stays at F_ = (paramagnetic phase). After the 
transition, it starts to roll down the potential hill in some direction 
[0(2) symmetry] which is set by the initial quantum (or thermal) 
fluctuations. 



responds to the re-heating process mentioned in Sec. Hi]) in a 
quantum analogue of the Kibble-Zurek ffl [l5ll mechanism, 
see also 0H1. 

Since the 0(2) invariance of the Hamiltonian ( fT3l ) with re- 
spect to rotations in the x,y-plane is broken by the ferromag- 
netic state, this phase supports topological (point) defects in 
the form of spin vortices in two spatial dimensions [29]. Such 
a defect corresponds to a spin configuration which cannot be 
deformed to a constant magnetization in a smooth way: Go- 
ing around a vortex in a circle in real space, the magnetization 
describes a circle in the internal space of the ferromagnetic 
ground-state manifold, i.e., the magnetization always points 
towards the defect (vortex) or away from it (anti-vortex). The 
number of vortices (minus the number of anti- vortices) within 
a given area is the winding number 9T. Calculating this wind- 
ing number for a circle of radius R, we obtain (for large R) the 
scaling law |12ll 



(m 2 ) ~R InR, 



(14) 



which lies in between the random phase walk model (predict- 
ing (Or) ~ R) and the random vortex gas model (predicting 
(9r) ~ R 2 ) frequently discussed in the literature. 



VIII. FINITE-SIZE SCALING 

So far, our studies were devoted to the thermodynamic limit 
of infinitely large systems (cf. the footnote in Sec. [V}. Strictly 
speaking, phase transitions are well-defined in this limit only. 
For systems of finite size, the sharp transition at a precisely lo- 
calizable critical point is typically broadened. Similarly, exact 
level crossings in the thermodynamic limit (as in Fig. [I} usu- 
ally become avoided level crossings for finite systems since 
boundary terms and finite-size effects lift the degeneracy in 
general. As a result, the response time may not diverge at the 
critical point but remain bounded (even though very large). In 
order to study non-equilibrium phenomena, the scaling of the 
response time (or the energy gaps) with system size is very 
important. 

As it turns out, this scaling may crucially depend on the 
order of the phase transition: For a typical first order transi- 
tion (cf. Fig.|4|, we have two competing ground states \\jf<) 



and | y>) which are locally distinguishable and become en- 
ergetically favorable just before \y<) = |Vb(<? < gent)} and 
after |tyA>) = \y/o(g > gait)) the transition. Due to their lo- 
cal distinguishability, their overlap decreases exponentially 
(y<| V>) ~ exp{— i^(n)} with system size n (e.g., number of 
lattice sites or spins). In the low-energy sub-space spanned 
by \Y<) and \Y>), we have four relevant matrix elements 
of the Hamiltonian. At the critical point g = g cnt , the diag- 
onal elements coincide (y<|i/|y<) = (y>|i/|y>), but the 
off-diagonal matrix elements (y < |/i r |y > ) = ((y>|H|y<))* 
lead to an avoided level crossing with the fundamental gap 
E\ —Eq = | (y<|J¥|y>) |. For a local Hamiltonian containing 
interactions of a limited number of sites (e.g., nearest neigh- 
bors) only H ~ poly(«), the off-diagonal matrix elements 
scale exponentially (y<\H \\jf>) ~ exp{— ^(n)} and hence 
the same scaling applies to the fundamental gap E\ — Eq ~ 
exp{— G{n)\. This result nicely matches the intuitive picture 
based on the energy landscape in Fig. |4] In order to stay in 
the ground state, the system has to tunnel through the energy 
barrier at the critical point. Since one would expect the height 
of the potential barrier to scale linearly with system size, the 
tunneling (i.e., response) time ~ 1 j(E\ — Eq) should increase 
exponentially lfl7ll . 

For transitions of second or higher order, on the other hand, 
such a tunnelling barrier is absent and there are no locally dis- 
tinct competing ground states at the critical point. Hence there 
is no reason for an exponential scaling of the response time in 
this situation. Indeed, for several analytically solvable exam- 
ples such as the quantum Ising model in a transverse field, the 
response time is found to scale polynomially (instead of expo- 
nentially) with system size lT30ll . Furthermore, the ground state 
changes less abruptly in second-order transitions than in those 
of first order, which also suggests that it should be easier to 
stay in the ground state. 



IX. ADIABATIC QUANTUM ALGORITHMS 

The intuition developed by considering the physics exam- 
ples above might also be useful for gaining a deeper under- 
standing in a rather different subject - quantum algorithms. 
Since the pioneering work of Shor and others, see, e.g., 
EH fl9tl . it is known that quantum computers (which oper- 
ate in the full Hilbert space of all quantum states) should be 
able to solve certain problems much faster than all (known) 
classical algorithms. However, actually constructing a quan- 
tum computer of sufficient size is extremely hard since the 
fragile superpositions necessary for operating the quantum al- 
gorithm are very vulnerable to decoherence caused by the in- 
evitable coupling to the environment. There are several ideas 
to overcome this problem: quantum error correcting codes, 
measurement-based ("one-way") quantum computers, noise- 
resistant (e.g., topological) qubit realizations etc. In the fol- 
lowing, we shall consider an alternative idea: adiabatic quan- 
tum algorithms 12011 . In this scheme, the solution to a given 
problem is encoded into the ground state of an appropriately 
constructed problem Hamiltonian H out . In order to arrive at 
the desired ground state of H out , one starts with an initial 
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FIG. 10: Interaction topology of the Hamitonian H6\ where each 
triple (a, j3 , y) corresponds to a triangle whose bonds govern the anti- 
ferromagnetic interaction M a a. 

Hamiltonian H m , whose ground state is known and can be pre- 
pared easily, and then slowly transform it to the final problem 
Hamiltonian H out , e.g., via a linear interpolation 

H(t) = [l-g(t)]H in + g(t)H out . (15) 

If this interpolation is slow enough, we end up in (or near) the 
final ground state encoding the solution to our problem 12111 . 
The computational complexity manifests itself in the runtime 
T, which is bounded from below by the response time deter- 
mined by the minimum gap. Now, for non-trivial problem 
Hamiltonians H out , there is typically a critical point where the 
fundamental gap E\ — Eq between the ground state and the 
first excited state becomes very small - which bears strong 
similarities to a quantum phase transition, see also 12211 . 

Let us study an explicit example. A rather interesting class 
of problems is exact cover-3, which corresponds to the fol- 
lowing requirement: Given a set of triples (a,/3,y) e N 3 , find 
at least one bit-string Za £ {0, 1} which satisfies the constraint 
Za +zp +Zy = 1 for all triples (a,j3,y). It can be shown that 
other tasks such as factoring can be mapped onto this problem 
class. A suitable problem Hamiltonian is given by the sum of 
(ct„ + Gp + + 1 ) 2 over all these triples, which yields 

n n 

Hout = M a,P°a (y p+ 2 L Na^a+const- (16) 

a,p = l a=l 

This corresponds to a frustrated anti-ferromagnet in an exter- 
nal field N a , where the interaction topology M a p is deter- 
mined by the set of triples, cf. Fig.fTOl 

Now, given the final problem Hamiltonian H oat above, what 
could be a suitable initial Hamiltonian //;„? Obviously, H m 
should not be diagonal in the rj^-basis - in this case, we would 
have exact (instead of avoided) level crossings and an adia- 
batic evolution would be impossible. Therefore, most previ- 
ous studies (see, e.g., [20]) adopted the form 

H m = Y L a< ■ (17) 

a 

On the other hand, inspired by the findings in Sec. IVIIII one 
could try to choose H m in a way such that the change form H m 
to //out is similar to a second-order phase transition. The ob- 
servation that symmetry-breaking/restoring phase transitions 
are often (though not always) of second order indicates one 
idea for doing this. An apparent symmetry of H out is its 0(2) 
invariance under rotations around the z-axis generated by 

J?=£o£. (is) 

a=\ 
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FIG. 1 1 : Plot of the runtime necessary for an adiabatic evolution for 
the conventional scheme d!7t and the XY network in Eqs. d!6t and 
J19l >. The hollow symbols correspond to especially hard problems 
within the class exact cover-3. 



cay and dephasing errors. However, this does not imply that 
these algorithms are completely free of decoherence: A finite 
temperature environment may heat up the system and cause 
excitations, which would result in computational errors (since 
the ground state, not an excited state, encodes the solution to 
our problem). Furthermore, non-equilibrium phenomena may 
also cause excitations at zero temperature due to the coupling 
to an environment. Since the energy gap typically becomes 
very small at some point of the interpolation, the algorithm 
will be very vulnerable near this point. Of course, the same 
applies to the critical point in more gen eral quantum phase 
transitions. For an adiabatic version ll23ll of the Grover lfl9ll 
quantum search algorithm - which corresponds to a first-order 
transition - we found that the induced excitation probability 
strongly depends on the spectral function /bath(w) of the en- 
vironment and scales as 112411 



/bath(A£ m i n ) 
A£Vnjn 



(20) 



Now, if we start with the ferromagnetic Hamiltonian 

H in = - £ M^o^ + o^), (19) 

a,j3 = l 

which is also 0(2) invariant, but whose degenerated ground 
states break this symmetry, one would expect that the resulting 
adiabatic quantum algorithm is very similar to a second-order 
phase transition, where the system finds it easier to stay in 
the ground state. In this way, physics intuition motivates a 
new way of constructing quantum algorithms lfl7ll . Indeed, 
numerical simulations suggest that the algorithm described in 
Eqs. ([ToT i and ( fT9l ) is superior (in average) to the conventional 
scheme ( fTTb , cf. Fig.fTTI 

X. DECOHERENCE 

One of the (original) motivations for adiabatic quantum 
computing was the robustness of the ground state against de- 



where AE m i„ is the minimum energy gap. Hence, if /b a th(<o) 
vanishes fast enough in the infrared (0 — > 0, the error rate can 
be kept under control. However, repeating the same analysis 
for the Ising model 112511 and the sweep through the second- 
order phase transition, we find that the excitation probability 
increases with system size n, independently of /bath (<*>)■ This 
suggests that second-order transitions are more suitable for 
adiabatic quantum computing - but, at the same time, more 
vulnerable to decoherence (which requires some error correc- 
tion techniques). 
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